Mobility Edge in Aperiodic Kronig-Penney Potentials with Correlated Disorder: 

Perturbative Approach 
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It is shown that a non-periodic Kronig-Penney model exhibits mobility edges if the positions of 
the scatterers are correlated at long distances. An analytical expression for the energy-dependent 
localization length is derived for weak disorder in terms of the real-space correlators defining the 
structural disorder in these systems. We also present an algorithm to construct a non-periodic but 
correlated sequence exhibiting desired mobility edges. This result could be used to construct window 
filters in electronic, acoustic, or photonic non-periodic structures. 
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The Kronig-Penney model has been widely used to ex- 
plore the characteristics of electrons in a periodic poten- 
tial, as this model provides one with perhaps the sim- 
plest instance of Bloch states. This model is also used 
systematically to provide estimates of the bandwidths in 
semiconductor superlattices with high reliability.o The 
Kronig-Penney model and its relation with superlat- 
tices has also been used in recent times to provide an 
implementation of the physics of random and quasi- 
periodic systemSjD and interesting experiments have been 
reported on arrangements such as the well-known Fi- 
bonacci sequences. u It is because of its importance and 
wide applicability that we focus our attention on the 
Kronig-Penney model. We will demonstrate that the 
aperiodic model with constant scattering potential but 
random spacings yields mobility edges if the disorder has 
long-range correlations, in sharp contrast to the situation 
for white noise potentials. 

The Kronig-Penney model in this study is given by a 
ID chain of delta-function scatterers with amplitude Un 
and centered at points z„. The Schrodinger equation for 
a particle moving in this random potential has the form, 

oc 

)■ (1) 



This is equivalent to the discrete equation 
[sin(/i. 



sin/Xn^iV'n+i -t-sin^„Vn-i = 
-I- ^J.n-l) + (Un/q) sin^„_i sin/x„] ?/'„, 



(2) 



where ipn = 4'{zn), fJ-n = q{zn+i ~ Zn), Q = Ve , and 
the energy is measured in units where h^/2'm = 1. The 
linear relation (^) between ipn-i, and V'n+i, can be 
easily arrived at by integrating Eq. (|l|) in the vicinity of 
sites n — 1, n, and 1, and substituting the amplitudes 
for the various constants in the piece-wise zero-potential 
regions between the scatterers. If the site potential is 
different from a delta function, a linear relation similar 
to Eq. (||) can be obtained using a general method de- 
scribed in Ref. |[ Therefore, Eq. (||) can be considered 



as a rather generic relation for ID chains with potential 
scatterers. We focus here on a sequence of scatterers of 
equal amplitude, C/„ = U, but with varying positions, 
i.e. the case of 'structural' or 'positional' disorder. An 
obvious experimental realization of this model is a semi- 
conductor superlattice with fluctuating period. 

The case of compositional disorder, i.e., a periodic ar- 
rangement of sites, z„ = n, with random amplitude [/„, 
has been studied intensively during the last decade. The 
importance of short-range correlations was first explored 
recentlyjj and a 'jsandom dimer' model was studied as a 
specific example.Ei Using this model, it was shown that 
short-range correlations in the infinite random sequence 
{Un} give rise to a discrete number of delocalized states, 
as well as to some anomalies in transport properties.Q 
The presence of such anomalies has been recently ob- 
served in experiments with GaAs-AlGaAs random-dimer 
superlattices.B Moreover, the non-trivial role of correla- 
tions in the formation of mobility edges has been pointed 
out by the study of localization in pseudorandom and 
incommensurate potentials.cl In contrast, only relatively 
recently the role of long-range correlations in random po- 
tentials has received special attention. The interplay of 
long-range correlations and disorder has been shown to 
lead to the existence of a continuum of extended states 
in the ene¥^£-|Spectrum and to the appearance of mobil- 
ity edges.E3llil These edges have been shown to exist in 
experiments of microwave transmission in a single-moda 
waveguide with a random array of correlated scatterers 
and their possible relevance for metal-insulator transi^ 
tions in 2D electron systems has been recently explored.113 

Localization length for weak disorder. To calculate the 
localization length 1{E) in the case^of structural disorder 
we use a Hamiltonian approacho In this scheme, the 
discrete Schrodinger equation is replaced by a classical 
Hamiltonian map for coordinate a;„ and momentum p„. 
For the case of Eq. (|^), conjugate variables are introduced 
as, Xn = ■>pn and p„ = (a;„cos^„_i - a;„_i)/ sin/i„_i. 
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Correspondingly, the discrete-time evolution of Xn and 
Pn is obtained from 



Pn+l = {Pn + An Xn) cos ^in " a;„sin/i„, 
Xn+1 ^ {Pn + A„ Xn) sin fj.„ + X„ COS ^„ . 



(3) 



This map describes the behavior of a linear rotator 
subjected to non-periodic delta-kicks with amplitude 
An = Un/q- Free rotation between kicks corresponds 
to free propagation between scatterers and each kick cor- 
responds to scattering at each 5-function potential. It is 
easy to check that the first equation in (0) is equivalent 
to ^ , while the second is reduced to an identity after p„ 
substitution. 

In what follows we consider the case of weak disorder, 
assuming that the deviation of the scatterers from their 
positions in a periodic lattice is small, |(5„| = q|z„ — n| <C 
1. We can then expand trigonometric functions in (j^) 
and up to second order in A„ — q{Sn+i — Sn), obtain the 
approximate map for constant kick amplitude A = U/q, 



Pn+l = cosq 
— sing 
Xn+l = cosq 
-I- sin (7 



A A 
Pri(l - -f ) + Xn{A - Xn - A^) 



PnX a + Ax 



PnXn + X„(l -I- AXn ~ 

A^ A^ 
P„(l - + Xn{A - Xn - A^) 



(4) 



this expansion is valid 



Note that since A„ cx g = 
only for low energies. 

In order to extract the effect that comes only from the 
positional non-periodicity, it is convenient to eliminate 
the mean field associated with the constant amplitude 
U . This can be done by a canonical transformation of 
variables, 



Pn 
Xn 



COS (j) Pn — a sin X„ 
sin (b Pn + a cos 6 Xn 



(5) 



The parameters of this transformation (a and (p) are ob- 
tained from the condition that to zeroth order in A„ 
the dynamical map for (P„,X„) be a simple rotation 



without kicks, i.e., P 



(0) 



= Pn^ cos 7 — Xn' sin 7, and 

-'^^+1 = Pn'^ sin 7 -I- Xn'^ cos 7. Applying this condition 
to Eqs. (|^) and (|j) we get after some algebra: 



(0), 



ri+l 



2(7, COS7 = cosg + (yl/2g) sing , . . 

-4 - 1 i 2^/[2sin9- A(l+cos9)]. ^ ' 



It is clear that 7 plays the role ofpthc Bloch number in 
the periodic Kronig-Penney model.u We can now rewrite 
the map (^) in terms of variables (P„,X„) and angle 7, 

Pn+l = (1 - A2 /2) {Pn cos 7 - X„ siu 7) 

-A„Q;2(P„sin7 + X„cos7), , , 

Xn+i - (1 - Xl/2) [Pn sin 7 + X„ cos 7) 
+A„a"^ [Pn cos 7 - X„ sin 7). 

Following our previous approach,0 it is convenient to 
introduce action-angle variables for map (0), 



Rn sin On, Xn = Rn COS 0„ 



(8) 



The inverse localization length for the original quantum 
model Eq. (|l|) (or equivalently, the Lyapunov exponent 
for the dynamical map in (^)), can be expressed in terms 
of the ratio Rn+i/Rn (see details in Ref. [l4|). 



l-\E) 



In 



n+l 



= In 



X 



ri+l 



X„ 



2 \ i?2 



(9) 



Here, angular brackets stand for the average over sites 
(kicks), so that, (...) = limAr^oo j;Y.n=i (■■■)■ Using the 
map (M) we can calculate the ratio Rn+i/Rn as 



(Rn+i/Rn) = 1 - f/A„ sin[2 (6'„ - 7)] /((7sin7) , . 
^2 r^-4sin2(^?„-7) + a4eos2(0„-7)-ll ^^^^ 



The fact that the ratio Rn+i/Rn is close to unity for 
small Xn is what motivated switching from the variables 
{pn,Xn) to {Pn,Xn) using (|^) . Thc logarithm in (^) can 
be expanded as ln(l + x) ~ x — x'^ /2, and, up to second 
order in A„, the average is performed over the unper- 
turbed motion given by the {Pi"\x^°^) variables, bmce 
this motion is a free rotation (in the old variables it is 
a rotation with periodic kicks), the angle variable 6n is 
clearly distributed uniformly within the interval [0,27r]. 
One then obtains that (^sin^ On) = {cos^ On) = 1/2, and 



l-\E) 



{^l)U' 
8q^ sin^ 7 



U 



qsm^ 



{Xn sm[2{dn - 7)]> 



(11) 



Thc first term in (|l^) gives the inverse localization length 
in an uncorrelated random potential. In this Born ap- 
proximation, it is proportional to the variance (A^) and 
to thc squared amplitude of the scattering potential, U^. 
Since (A^) cx g^, the factor q^ disappears and the only 
energy dependence is due to the factor sin^ 7 in the de- 
nominator. At the edges of the allowed zones sin 7 — 0, 
and here the localization length 1{E) approaches zero. 
A similar enhancement of localization in the vicinity of 
the band edgeslijii3 has stimulated the study of photonic- 
band-gap materials in the last decade. 

The second term in ( |ll| ) describes the contribution 
of correlations in thc scattering potential. To cal- 
culate exphcitly the correlator (A„ sin[2(0„ — 7)]), one 
needs the recursion relation for the angle variable 
On- Since this correlator already contains a factor A„, 
only linear terms in the recursion relation are needed 
from Eqs. (^ and (||), so that On = On-i — 7 — 
A„_i [a^ — U siv?{On^i ~ 7)/(gsin7)]. The correlator 
(A„ sin[2(^„ — 7)]) can be written as a Fourier series in 
Bloch number 7, where the dipiensionless correlators S,{k) 
are the Fourier coefRcientsO 



(A„ sin[2(0„ - 7)]) = -t;^ ^('^) cos(27n) . (12) 
2 sm 7 
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Here = (A^), with ^k) - (A„+,.A„) /A^, and 
= ^n+i — ^n- Note that the locahzation length is de- 
termined by the statistical properties of the sequence of 
relative displacements A„, and not by the displacements 
dn themselves. 

Substituting the correlator (^2|) into (|ll]) we get the 
final result for the inverse localization length, 

1 [/^A^ °° 
Jr^ = „ . 2 fil) = 1 + 2 V cos(27n). 



n=l 



(13) 



This formula has the same structure as that obtained for 
a tight-binding (and the corresponding Kronig-Penney) 
model wi th . ra ndom amplitudes U„, but equidistant sites 
(z„ — 7i).EJii3 These three different models have a differ- 
ent dependence on energy via the factors and sin^ 7. 
The present case exhibits the weakest dependence of the 
localization length on energy within the allowed zone 
since a factor = E appears in other cases but not 
here. This property should be favorable for the experi- 
mental observation of mobility edges in superlattices with 
positional disorder. 

Correlations and mobility edge. If the sequence of ran- 
dom displacements A„ is uncorrelated, ^{k) = 0, the 
localization length is given by the first term in (p^. In 
the opposite limit of a completely correlated sequence, 
^(fc) =const., the displacements are independent on the 
site number, (5„ = Sq, giving a regular sequence with pe- 
riod l-l-(5o, and extended states, l^^ = 0, for all energies. 
A smooth transition between these two limits can be de- 
scribed by an exponential function, ^(fc) = exp(— fc/fco), 
where kg is a correlation radius. Substituting this form 
into ([l3|) one obtains 



l-\E) 



sinh (1/fco) 



8 sin 7 cosh(l/fco) - cos(27) 



(14) 



For any finite fco all states are localized. Only in a 
periodic lattice (fco = 00) does the inverse localization 
length ([1^ ) vanish and the states become delocalized. 
This localization-delocalization transition occurs simul- 
taneously for all energies, and a mobility edge does not 
appear in the spectrum. This conclusion is valid, for ar- 
bitrary amplitude of the scattering potential [/EO A dis- 
crete number of delocalized states can in fact appear if 
the binary correlatac,^(fc) oscillates with exponentially 
decaying amplitude.li3 This can be obtained from (^^, if 
the correlation radius fco is allowed to take complex val- 
ues. On the other hand, a mobility edge may appear if 
correlations-jdpcay not exponentially but according to a 
power law.E2rt3 We show numerically below that mobil- 
ity edges exist if, e.g., ^(fc) oc 1/fc. 

Designed mobility edges. In a real (or numerical) ex- 
periment one needs to know explicitly the displacements 
Sn which provide a desirable dependence 1{E). This leads 
us to the 'inverse problem' in the theory of localization. 



The solution would give us in general a random scat- 
tering potential {C/„,(5„} through the dependence 1{E). 
In the present case of positional disorder, [/„ = U, we 
need to calculate only the relative displacements. One 
can explicitly evaluate the correlator ^ (fc) as the Fourier 
coefficient, 



Tr/2 



ifiij) cos(2fc7)(i7 , 



(15) 



where the energy dependence of </'(7) = 
Ssin^ ^/{U^ 1(E)), is assumed to be known. The 
energy E ~ is expressed through 7 via the disper- 
sion relation in (||). It is easy to check that the binary 
correlator of a sequence (A„/A) = X]fc°=-oo /^(^) ^n+k , 
coincides with ^(fc) if are random numbers with zero 
mean an| d | pait | variance, and where the function /3(fc) is 
given bytiltSHla 

2 /"^/s 

/3(fc) = -/ v/R7)cos(2fc7)d7. (16) 

7^ Jo 

Once the relative displacements A„ are known, a se- 
quence of absolute displacements can be easily calcu- 
lated, setting So = 0, for example, and then di ~ Aq, 
62 = Aq-I-Ai, . . . , 6n = X]fc=o ^1^- This procedure allows 
the calculation of displacements S„ for any energy depen- 
dence of the localization length 1{E), including situations 
with mobility edges. Appropriately correlated elements 
may be used for fabrication of effective filters of electric 
cal or optical signals, even if the system is not periodic.t2l 
The bandwidth of a filter can be made arbitrarily wide 
or narrow depending on the statistical properties of the 
random sequence used (via the function (3{k)). 

Numerical examples. In order to examine our predic- 
tions, we construct explicit random sequences {i5n} for 
which the function 1^9(7) in (|l^) has four mobility edges 
at 7i = 0.27ri, i = 1, . . . , 4. The positions of the mobility 
edges are chosen within the interval < 7 < tt, symmet- 
rically about 7 — tt/2. Via relation (|^), the position of 
the mobihty edges on the axis q/tt are given by 0.326, 
0.478, 0.649, and 0.850, for the mean field amphtude 
U = 0.7. In the ideal lattice with the same strength of the 
potential, the first allowed band lies between qi/n « 0.26 
(7 = 0) and qr/'n « 1 (7 = tt). Numerical data are 
given for two complimentary situations ^pi^2{l)'- fi van- 
ishes (i.e. 1(E) — 00 and states are delocalized) in the 
region 7 S (0,7i) U (72,73) U (74, tt); while (p2 vanishes 
in the complementary region 76 (71, 72) U (73, 74). Out- 
side these regions, the function (/3i,2(7) is a constant de- 
fined by the normalization condition ^1,2(0) = 1. From 
Eq. (p^), we get that because of the presence of sharp 
mobility edges the correlators (Fourier components of a 
discontinuous function) decay slowly: 



a(fc) = -i.56(fc) = 

(5/27rfc) [sin(0.87rfc) - sin(0.47rfc)] 



(17) 
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FIG. 1. Inverse localization length for band with four mo- 
bility edges, three regions of extended states, and different 
amplitude of disorder: circles (o), A = 0.05; crosses (x), 
A = 0.15. Solid lines show Eq. @ for A = 0.05. Inset 
shows TV — 10^ sites averaged over only 100 samplings. 




q/71" 



FIG. 2. Results for complementary system to Fig. 1, with 
four mobility edges but only two regions of extended states. 
Lines through circles (A = 0.05) and crosses (A = 0.15) are 
guides to the eye. 

We show the corresponding data in Fig. 1 and 2 for 
U — 0.7. The analytical dependence (|l^) for the dimen- 
sionless inverse localization length A = 8/l{E)U'^A'^ is 
shown by the full lines for A = 0.05 in the figures. Nu- 
merical data are obtained for a large sample size N — 10^, 
with two amplitudes of disorder, A = 0.05 (circles) and 
A = 0.15 (crosses). The insets show results for a much 
shorter sample, N = 1000, with additional average over 
100 different realizations of disorder but the same corre- 
lations. One can see that for small disorder, A = 0.05, 
Eq.(|l3[) describes the numerical results very well, with 



only minor deviations close to the onset of mobility qi 
and band edges, qi^r ■ Notice that Eq.(|l^) is excellent 
at giving the mobility edges at the prescribed energies, 
which proves its usefulness, and breadth of applicability. 
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versity. 
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The Kronig-Penney model has been widely used to ex- 
plore the characteristics of electrons in a periodic poten- 
tial, as this model provides one with perhaps the sim- 
plest instance of Bloch states. This model is also used 
systematically to provide estimates of the bandwidths in 
semiconductor superlattices with high reliability.^ The 
Kronig-Penney model and its relation with superlat- 
tices has also been used in recent times to provide an 
implementation of the physics of random and quasi- 
periodic systems,^ and interesting experiments have been 
reported on arrangements such as the well-known Fi- 
bonacci sequences.^ It is because of its importance and 
wide applicability that we focus our attention on the 
Kronig-Penney model. We will demonstrate that the 
aperiodic model with constant scattering potential but 
random spacings yields mobility edges if the disorder has 
long-range correlations, in sharp contrast to the situation 
for white noise potentials. 

The Kronig-Penney model in this study is given by a 
ID chain of delta-function scatterers with amplitude Un 
and centered at points z„ . The Schrodinger equation for 
a particle moving in this random potential has the form, 

oo 

{n^/2m)ij"{z) + Eij{z)= J2 Uni^{z„)S{z - z„) . (1) 

n = — oo 

This is equivalent to the discrete equation, 

[sin{nn + Hn-i) + {Un/q) sin Hn-i sin Hn]lpri, 

where ip„ = ip{z„), Hn = q{zn+i - z„), q = ^/E, and 
the energy is measured in units where Ti^ /2m = 1. The 
linear relation (2) between ipn-i, i'n, and ipn+i, can be 
easily arrived at by integrating Eq. (1) in the vicinity of 
sites n — 1, n, and n-\-l, and substituting the amplitudes 
for the various constants in the piece-wise zero-potential 
regions between the scatterers. If the site potential is 
different from a delta function, a linear relation similar 
to Eq. (2) can be obtained using a general method de- 
scribed in Ref. 4. Therefore, Eq. (2) can be considered 



as a rather generic relation for ID chains with potential 
scatterers. We focus here on a sequence of scatterers of 
equal amplitude, Un = U, but with varying positions, 
i.e. the case of 'structural' or 'positional' disorder. An 
obvious experimental realization of this model is a semi- 
conductor superlattice with fluctuating period. 

The case of compositional disorder, i.e., a periodic ar- 
rangement of sites, z„ = n, with random amplitude Un, 
has been studied intensively during the last decade. The 
importance of short-range correlations was first explored 
recently,® and a 'random dimer' model was studied as a 
specific example.® Using this model, it was shown that 
short-range correlations in the infinite random sequence 
{Un} give rise to a discrete number of delocalized states 
as well as to some anomalies in transport properties.^ 
The presence of such anomalies has been recently ob- 
served in experiments with GaAs-AlGaAs random-dimer 
superlattices.® Moreover, the non-trivial role of correla- 
tions in the formation of mobility edges has been pointed 
out by the study of localization in pseudorandom and 
incommensurate potentials.^ In contrast, only relatively 
recently the role of long-range correlations in random po- 
tentials has received special attention. The interplay of 
long-range correlations and disorder has been shown to 
lead to the existence of a continuum of extended states 
in the energy spectrum and to the appearance of mobil- 
ity edges. "'^'^'"'^"'^ These edges have been shown to exist in 
experiments of microwave transmission in a single-mode 
waveguide with a random array of correlated scatterers, "'^^ 
and their possible relevance for metal-insulator transi- 
tions in 2D electron systems has been recently explored. 

Localization length for weak disorder. To calculate the 
localization length 1{E) in the case of structural disorder 
we use a Hamiltonian approach.^'* In this scheme, the 
discrete Schrodinger equation is replaced by a classical 
Hamiltonian map for coordinate Xn and momentum pn ■ 
For the case of Eq. (2), conjugate variables are introduced 
as, Xn = tpn and pn = [xncosnn-i " «n - 1 ) / sin _ i . 
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Correspondingly, the discrete-time evolution of x„ and 
Pn is obtained from 



p„+i = (p„ +A„x„)cosfi„ - x„smfi„, 
x„+i = {pn + A„ x„) sin iJ,„ + x„cosii„. 



(3) 



This map describes the behavior of a linear rotator 
subjected to non-penodtc delta-kicks with amplitude 
An = Un/g- Free rotation between kicks corresponds 
to free propagation between scatterers and each kick cor- 
responds to scattering at each (5-function potential. It is 
easy to check that the first equation in (3) is equivalent 
to (2), while the second is reduced to an identity after p„ 
substitution. 

In what follows we consider the case of weak disorder, 
assuming that the deviation of the scatterers from their 
positions in a periodic lattice is small, \S„ \ = q\z„ — n\ <^ 
1. We can then expand trigonometric functions in (3) 
and up to second order in X„ = q{S„-\.i — S„), obtain the 
approximate map for constant kick amplitude A = U/q, 



Pn + l 



COS q 
■ sin q 



Pn(l - ^) + Xn{A- Xn - A^\ 



PnK + «n(l + AXr^ 



[-n + l 



cosg ^PnX„ + x„(l + AX„ f) 

+ sin q \p„(l - ^) + x„(A - X„ - A^) 



(4) 



Note that since Xn <x q = V E , this expansion is valid 
only for low energies. 

In order to extract the effect that comes only from the 
positional non-periodicity, it is convenient to eliminate 
the mean field associated with the constant amplitude 
U. This can be done by a canonical transformation of 
variables, 



p„ = a 1 cos (j> P„ — a sin (j> X„ 
x„ = sin (j> P„ + a cos (j> X„ 



(5) 



The parameters of this transformation (a and (f)) are ob- 
tained from the condition that to zeroth order in X„ 
the dynamical map for {P„,X„) be a simple rotation 
without kicks, i.e., -P^!|_\ = -Pi'^'' cos 7 — X^'^'' sin 7, and 

x'f^-^ = P^f ^ sin 7 -|- xi^"^ cos 7. Applying this condition 
to Eqs. (5) and (4) we get after some algebra: 

<f> = ^q , cos 7 = cos q + [A/2q) sin q . 



l + 2A/\2sinq - A{1 + cos q) 



(6) 



It is clear that 7 plays the role of the Bloch number in 
the periodic Kronig-Penney model. ^ We can now rewrite 
the map (4) in terms of variables {P„,X„) and angle 7, 



P„+i = (1 - Xl/2) (P„ cos 7 - X„ sin 7) 

-X„a^{Pn sin 7 -|- X„ cos 7), 
X„+i = {I -Xl/2) (P„sin7-FX„cos7) 
+X„a~^{P„ COS7 - X„ sin 7). 



(7) 



Following our previous approach, it is convenient to 
introduce action-angle variables for map (7), 



P„ = Rn sin t 



X„ = Rn COS ( 



(8) 



The inverse localization length for the original quantum 
model Eq. (1) (or equivalently, the Lyapunov exponent 
for the dynamical map in (3)), can be expressed in terms 
of the ratio R„^i/R„ (see details in Ref. 14), 



In 



n + l 



In. 



X, 



n + l 



X„ 



1 / Rl+i 



(9) 



Here, angular brackets stand for the average over sites 
(kicks), so that, (...) = limAr^oo -kYln=i (•••)• Using the 
map (7) we can calculate the ratio R„+i/Rn as 



[Rn+i/Rnf = l-UXn sin [2 {e^ - 7)] /{q sin 7) 



a ■* sin^ 



{On - 7) + ' 



7)-l] 



(10) 



The fact that the ratio R„^i/R„ is close to unity for 
small X„ is what motivated switching from the variables 
{Pn,Xn) to {P„,X„) using (5). The logarithm in (9) can 
be expanded as ln(l + x) x — x^ /2, and, up to second 
order in X„ , the average is performed over the unper- 
turbed motion given by the {P!l°\xi°^) variables, bmce 
this motion is a free rotation (in the old variables it is 
a rotation with periodic kicks), the angle variable 6„ is 
clearly distributed uniformly within the interval [0,27r]. 
One then obtains that (sin^ On) = (cos^ 9n) = 1/2, and 



U 



8g2sin^7 gsin7 



{X„sin[2{0„-j)]) . (11) 



The first term in (11) gives the inverse localization length 
in an uncorrelated random potential. In this Born ap- 
proximation, it is proportional to the variance {X'l) and 
to the squared amplitude of the scattering potential, . 
Since {X'^) cx g^, the factor disappears and the only 
energy dependence is due to the factor sin^ 7 in the de- 
nominator. At the edges of the allowed zones sin 7 = 0, 
and here the localization length 1{E) approaches zero. 
A similar enhancement of localization in the vicinity of 
the band edges^®'^® has stimulated the study of photonic- 
band-gap materials in the last decade. 

The second term in (11) describes the contribution 
of correlations in the scattering potential. To cal- 
culate explicitly the correlator (A„ sin[2(6'„ — 7)]), one 
needs the recursion relation for the angle variable 
6n ■ Since this correlator already contains a factor X„ , 
only linear terms in the recursion relation are needed 
from Eqs. (7) and (8), so that 6„ = 0„-i — 7 — 
A„_i [a^ — [/sin^(6'„_i — 7)/(g sin7)] . The correlator 
(A„ sin[2(6'„ — 7)]) can be written as a Fourier series in 
Bloch number 7, where the dimensionless correlators £,{k) 
are the Fourier coefficients,^^ 



(A„ sin[2{0„ - j)]) 



UqA^ 
2 sin 7 



^e(n)cos(27n). (12) 
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Here = (A^), with ^{k) = (A„+feA„)/A2, and 
A„ = S„+i — S„. Note that the localization length is de- 
termined by the statistical properties of the sequence of 
relative displacements A„, and not by the displacements 
S„ themselves. 

Substituting the correlator (12) into (11) we get the 
final result for the inverse localization length, 

1 IPA^ °° 
TTin = Q ■ 2 ^(7) = l + 2^^(n)cos(27n). 

l(-tj) 8sm 7 

^ ' ' n=l 

(13) 

This formula has the same structure as that obtained for 
a tight-binding (and the corresponding Kronig-Penney) 
model with random amplitudes Un , but equidistant sites 
(z„ = n)}^'^'^ These three different models have a differ- 
ent dependence on energy via the factors and sin^ 7. 
The present case exhibits the weakest dependence of the 
localization length on energy within the allowed zone 
since a factor = E appears in other cases but not 
here. This property should be favorable for the experi- 
mental observation of mobility edges in superlattices with 
positional disorder. 

Correlations and mobility edge. If the sequence of ran- 
dom displacements A„ is uncorrelated, ^(fc) = 0, the 
localization length is given by the first term in (13). In 
the opposite limit of a completely correlated sequence, 
^(k) =const., the displacements are independent on the 
site number, S„ = Sq, giving a regular sequence with pe- 
riod I + Sq, and extended states, l~^ = 0, for all energies. 
A smooth transition between these two limits can be de- 
scribed by an exponential function, ^(fc) = exp(— fc/fco), 
where kg is a correlation radius. Substituting this form 
into (13) one obtains 

smh(l/fco) 

^ ' 8sin2 7 cosh(l/fco) -cos(27) • ^ ' 

For any finite kg all states are localized. Only in a 
periodic lattice {kg = 00) does the inverse localization 
length (14) vanish and the states become delocalized. 
This localization-delocalization transition occurs simul- 
taneously for all energies, and a mobility edge does not 
appear in the spectrum. This conclusion is valid for ar- 
bitrary amplitude of the scattering potential U .^"^ A dis- 
crete number of delocalized states can in fact appear if 
the binary correlator ^(fc) oscillates with exponentially 
decaying amplitude.^® This can be obtained from (14), if 
the correlation radius kg is allowed to take complex val- 
ues. On the other hand, a mobility edge may appear if 
correlations decay not exponentially but according to a 
power law.^*^^^^ We show numerically below that mobil- 
ity edges exist if, e.g., ^(fc) cx l/k. 

Designed mobility edges. In a real (or numerical) ex- 
periment one needs to know explicitly the displacements 
S„ which provide a desirable dependence 1{E). This leads 
us to the 'inverse problem' in the theory of localization. 



The solution would give us in general a random scat- 
tering potential {Un,S„} through the dependence 1{E). 
In the present case of positional disorder, (/„ = U , we 
need to calculate only the relative displacements. One 
can explicitly evaluate the correlator ^(fc) as the Fourier 
coefficient, 

m = - I ' ^(7)cos(2fc7)rf7, (15) 

where the energy dependence of = 
8 sin^ 7/([/^A^/(£')), is assumed to be known. The 
energy E = q'^ \s expressed through 7 via the disper- 
sion relation in (6). It is easy to check that the binary 
correlator of a sequence (A„/A) = ^^.^o Pi^) Zn+k , 
coincides with ^(fc) if Zk are random numbers with zero 
mean and unit variance, and where the function l3{k) is 
given by-'--'-'-'-^'-'-^ 

2 r/^ , 

m = - ^Ml) cos{2kj)dj . (16) 

^ Jo 

Once the relative displacements A„ are known, a se- 
quence of absolute displacements can be easily calcu- 
lated, setting So = 0, for example, and then Si = Aq, 
S2 = Aq-I-Ai, . . . , Sn = Yl2=o '^k- This procedure allows 
the calculation of displacements S„ for any energy depen- 
dence of the localization length 1{E), including situations 
with mobility edges. Appropriately correlated elements 
may be used for fabrication of effective filters of electri- 
cal or optical signals, even if the system is not periodic. 
The bandwidth of a filter can be made arbitrarily wide 
or narrow depending on the statistical properties of the 
random sequence used (via the function l3{k)). 

Numerical examples. In order to examine our predic- 
tions, we construct explicit random sequences {S„} for 
which the function ip{'j) in (13) has four mobility edges 
at 7j' = 0.27ri, i = 1, . . . , 4. The positions of the mobility 
edges are chosen within the interval < 7 < tt, symmet- 
rically about 7 = 7r/2. Via relation (6), the position of 
the mobility edges on the axis q/ir are given by 0.326, 
0.478, 0.649, and 0.850, for the mean field amplitude 
U = 0.7. In the ideal lattice with the same strength of the 
potential, the first allowed band lies between qi/n « 0.26 
('J = 0) and qr/TT ~ 1 (7 = tt). Numerical data are 
given for two complimentary situations ipi^2{l)'- f\ van- 
ishes (i.e. 1{E) = cxD and states are delocalized) in the 
region 7 G (0,7i) U (72,73) U (74, tt); while vanishes 
in the complementary region 7 G (71, 72) U (73, 74). Out- 
side these regions, the function f>\^'i{l) is a constant de- 
fined by the normalization condition ^1^2(0) = 1. From 
Eq. (15), we get that because of the presence of sharp 
mobility edges the correlators (Fourier components of a 
discontinuous function) decay slowly: 

6(fc) = -1.56(fc)= (17) 
(5/27rfc) [sin(0.87rfc) - sin(0.47rfc)] . 
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FIG. 1. Inverse localization length for band with four mo- 
bility edges, three regions of extended states, and different 
amplitude of disorder: circles (o), A = 0.05; crosses (x), 
A = O.fS. SoUd lines show Eq. (f3) for A = 0.05. Inset 
shows = fO^ sites averaged over only fOO samplings. 
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FIG. 2. Results for complementary system to Fig. f , with 
four mobility edges but only two regions of extended states. 
Lines through circles (A = 0.05) and crosses (A = 0.f5) are 
guides to the eye. 

We show the corresponding data in Fig. 1 and 2 for 
U = 0.7. The analytical dependence (13) for the dimen- 
sionless inverse localization length A = 8/l{E)U^ A"^ is 
shown by the full lines for A = 0.05 in the figures. Nu- 
merical data are obtained for a large sample size N = 10®, 
with two amplitudes of disorder, A = 0.05 (circles) and 
A = 0.15 (crosses). The insets show results for a much 
shorter sample, N = 1000, with additional average over 
100 different realizations of disorder but the same corre- 
lations. One can see that for small disorder, A = 0.05, 
Eq.(13) describes the numerical results very well, with 



only minor deviations close to the onset of mobility rji 
and band edges, qi r ■ Notice that Eq.(13) is excellent 
at giving the mobility edges at the prescribed energies, 
which proves its usefulness, and breadth of applicability. 
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